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Abstract 

^ ' A numerical search for straight superconducting vortices in a U{1) 

■ model with a Ginzburg-Landau potential containing a cubic term, is pre- 
^ , sented. Such vortices exist in a small numerically determined region. The 
^SJ I reasons of their existence in that narrow region of the parameter space, as 
O . well as of their instability in the rest of the parameter space, are explained. 
^ . Then, the results of a numerical search for axially symmetric solitons in a 

U{1) X U{1) model with higher derivative terms, which is based on |18| . 
p • are presented and discussed. 

■ 1 Introduction 



X 



The evolution of the Universe is believed to involve several symmetry breaking 



;h ' phase transitions, out of which, topological defects, such as strings, are created 

^ [HE]- These transitions can be examined in the framework of condensed matter 

systems. Although there are differences from the cosmological case where, rela- 
tivistic dynamics must be used and gravity is important, the formation of such 
defects in the laboratory [3|, can provide helpful hints for cosmology. Topolog- 
ically stable knots and vortex-like structures in general, are of wide interest in 
condensed matter physics. For example, one can think of Bose-Einstein conden- 
sates (BEC) (i.e. see [4]), vortices in superfluid Helium-3 and Helium-4 [5] , or 
nematic liquid crystals f6l[7]- 

Also, in the framework of high energy physics, future experiments in LHC 
could answer whether metastable particle-like solitons exist in minimal supersym- 
metric Standard Model or two-Higgs Standard Model (2HSM) or not. In [8]-[Tn]. 
work on classically stable, metastable quasi-topological domain walls and strings 
in simple topologically trivial models, as well as in the 2HSM has been done. 
These solutions are local minima of the energy functional and can quantum me- 
chanically tunnel to the vacuum, not being protected by an absolutely conserved 
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quantum number. One can also find other interesting subjects involving super- 
conducting vortex rings such as, rotating superconducting rings electroweak 
strings [T2l [T3] or work on such rings in SU (2) non-Abelian Yang-Mills-Higgs 
model [H]. Finally, twisted semilocal vortices examined in [Hj can be connected 
to the models we present below, while one can also search if stable rings can exist 
in that model. 

This paper consists of two parts. In the first part, we consider a U{1)a model 
with a modified Ginzburg-Landau (GL) potential. The modification has to do 
with the addition of a cubic term. In Thermal Field Theory, such term comes 
from the 1-loop radiative corrections to the GL potential [I6|. We search whether 
this model can admit stable strings or not. The features of such strings, if they 
exist, are the supercurrent which fiows on the surface of the defect within a 
certain finite width, as well as, a magnetic fiux in the interior of the defect. This 
magnetic fiux is a consequence of the existence of the supercurrent. One has 
to keep in mind that the magnetic field can penetrate in a certain depth inside 
the superconducting regions where the supercurrent fiows. If the penetration 
depth is greater than the width of the superconducting surface, then the defect 
becomes unstable and can be destroyed. The GL potential we use here, is used in 
condensed matter physics as well (see [T7j and references therein) thus, one can 
also make interesting connections with that sector and the experiments described 
in the first paragraph. 

Stable defects of this U{1)a model, can also be used to form torus-like strings 
and study their stability. This can happen by taking a piece of such straight string 
and periodically connect its ends together. These string loops are examined in 
[T8] but in the framework of a U{1) x U{1) model, where the existence of the 
defect is ensured for topological reasons [HIES]. That model, is a continuation 
of previous work [8]-[l0]. In the second part of this paper, we examine a modified 
version of the model in [1^. We add higher derivative terms which might help in 
stabilizing the torus-shaped soliton. We investigate that possibility, present and 
analyse our results. 

2 The U{1)a model 

This model consists of a complex scalar field ip and a gauge field A^. The La- 
grangian density describing our system is: 



where the covariant derivative is = + ieA^, the strength of the field is 
F'^iu = df^Ai, — dyA^, while e is the U{1)a charge. We choose the potential 
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Figure 1: The relative position of the supercurrent as well as of the magnetic field on 
a j;y-profile of the system on the right. The left picture is strongly reminiscent of an 
infinite solenoid. 



where a, /3,7 constants. We can set 7 = 1. Since has dimensions of mass, 
we count the energy of the system in units of ^7. The vacuum is \ip\ = 0. This 
vacuum leaves unbroken the gauge symmetry U{1)a- When 



m = m\ = ^ (3) 

(for P > 2), we have U{1)a 1 giving non-zero mass to A. Thus, one may 
generate an electric current flowing along regions where |^| 7^ 0. In figE] one 
can see the shape of the potential. The eq.(l3|) gives the position of the minimum 
of interest for every f5 > 2. When (3 = 2, the secondary (non-trivial) minimum 
of the potential disappears, at j-^l = 1 position. When 2 < (3 < = /Scrit, it 
becomes zero only at \i(j\ = 0, while another minimum with non-zero forms. 
When P = Pcritj the potential has another zero at \ip\ = \/2 which is also a local 
minimum. Finally, when /? > Pcrit it has two more zeros at 



between which, it becomes negative (figEj). The mass spectrum is 



a 



2.1 The U{1)a model: Search for stable vortices 

We are interested in conflgurations having cylindrical symmetry, that is, infinite 
straight strings (figH]) . The field ip can be non- vanishing on a cylindrical surface 
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of specific radius. At infinity (p — > oo), we have the vacuum of the theory = 0. 
The ansatz for the fields is: 



(6) 



where M the winding number of the field ijj and p, z are the cylindrical 
unit vectors. We use cylindrical coordinates {t,p,(p,z), with space-time metric 
9ij.i^ — diag{l, —1, —p^, —1). We work in the = gauge. For the gauge field we 
suppose the above form based on the following thought: The A field is the one 
produced by the supercurrent fiowing on the cylindrical surface. The current is 
in the (p direction thus, we expect the non- vanishing component to be A^ and the 
amplitude P of ^ to be independent of ip. As it concerns the scalar field ip, since 
it follows the geometry of the cylindrical defect, we expect that its amplitude is 
independent of z as well. 

With the above ansatz, the energy functional for minimization takes the form 



E^2n pdp 



7^,{d,A^r + {d,P? + ^{eA^ + Mf + U{P) 
z,p p 



and the potential is 



(7) 



(8) 



The gauge field A has a magnetic field of the form 



p dp 



while, the field equations are 



dlP + -dpP - -{eA^ + Mf -'^(p^-(3P + l^=Q 

dlA^ - -dpA^ - 2eP\eA^ + M) = 
The usual rescaling arguments lead to the virial relation: 





h = 27r 




-27r 



pdp—— — — + - 



(9) 

(10) 
(11) 

(12) 
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Figure 2: The potential for a = 1 and /3 = 2, /? = Pcrit = (dashed line) and (3 = 2.17 
(bottom line). The plot is U vs. P. 



For a solution of the model, we theoretically must have Ii + I2 = 0. In fact, we 
define the index V = . We want this index as small as possible. Other 

virial relations can be found as follows. For example, one can consider the double 
rescaling of p — > Ap and either P — > fiP or fiA^ or even both of the fields 

and then demand dxE\x=i=ij, = = d^E\x=i=f_i- 



2.2 The U{1)a model: Numerical results 

We use a standard minimization algorithm to minimize the energy functional (I7j) . 
The algorithm is written in C. One can find details about the algorithm used, on 
page 425 of [22] but, briefiy, the basic idea is this: Given an appropriate initial 
guess, there are several corrections to it, having as a criterion the minimization 
of the energy in every step. When the corrections at the value of the energy 
are smaller than ^ 10~^ the program stops and we get the final results. We 
are interested in final configurations having non-trivial energy. This signals the 
existence of a stable vortex with that energy. We check our results through virial 
relation (fT2l) . Finally, our results must also satisfy the field equations (fTOl) . 
The initial guess we use for our computation is: 

Pip) = ^ip"'(l-tanh(0.2p2)) 
A^ = -tanh(^2P^) 

where ^1, .^2 are constants, the value of which, depends also on the location of 
the minimum (say = \ipo\) of the potential. For the final configurations we 
present on the figures, we chose ^1 = 1.75, ^2 = 2. 
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Figure 3: In the first two graphs we have the initial guess (dashed hues ) as well 

as the final configuration of fields (solid lines — ) P{p), Aip{p). We chose M = 1, 
e = 1, P = 2.17, a = 43.7. The energy E = 35.3 and virial is 10~^. The bottom graph 
gathers all the fields. For the field P, the area between the horizontal lines P3 and P4 
is energetically favorable (see also figE]). 

The initial guess also satisfies the appropriate asymptotics 

• near p = 0: P ~ p^' , ~ 

• at infinity: P ^ 0, exponentially 

while P must be non-zero somewhere between p = and p ^ 00. One must 
make a careful choice of the initial guess. That is to say, the maximum value of 
P at the initial guess (dashed lines) must be inside the favorable area denoted by 
the horizontal lines in figslSH This area is dictated by the form of the potential 
and especially by its negative sectors (see figE]). 

Since, for values of P where 2 < /? < Pcru is valid, we find no non-trivial 
solution, we searched and tried to find out what happens when we make the 
minimum of interest deeper. For that reason we searched in the region where 
P > Peril- It is possible to find solutions to this model until [3 reaches 2.13 (from 
above). Under this value, this is difficult if not impossible. Even at /3 = 2.13 
(figS]) we use great values of a in order to find the solution exhibited. The 
relation between a and [3 can be found in figEl The solutions of the model for 
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Figure 4: In the first two graphs we have the initial guess (dashed hues ) as well as 

the final configuration of fields (solid lines — ) P{p), A^{p). We chose M = 1, e = 1, 
/5 = 2.13, a = 1104. The energy E = 115.2 and virial is 2 • 10"^. The bottom graph 
gathers all the fields. For the field P, the area between the horizontal lines Pi and P2 
is energetically favorable (see also figE]). 

two different values of /? > Pcrit, are shown in figslMlland a comparison betv^^een 
them in figl5] in order to observe their different features. 

2.2.1 Analysis for 2.13 </? < 2.2 

We observe that for a specific /3, there is a small range for the parameter a of the 
potential, where the model exhibits the solution presented in the figures. Out of 
this small range and for greater a, we end up to a final configuration of negative 
energy. This happens because, the bigger the parameter a becomes, the stronger 
the potential is, thus the field P strongly prefers to acquire the value where the 
non-trivial minimum of the potential is (see figE]), in order to decrease further 
the energy. But, for /? > Pcru we have ?7 < at the position of that minimum, 
which also enforces the total energy E to be negative in this case. We have to 
note that, virial relation of such a final configuration is not satisfied due to the 
great values the term {dpPy acquires around p = and p ^ 00. This is clear if 
one directly observe (fT2l) . which can not be satisfied for ?7 < 0. 

On the other hand, for smaller values of a, we get the trivial configuration 
(P = 0, = 0), as the benefit from the potential term is no longer satisfactory 
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Figure 5: Here we have a plot of the above sets of parameters in order to compare 
the changes of P{p) and Ba{p) as a goes from 2.17 (dashed and dotted - • - ), to 2.15 
(dashed ) and 2.13 (dotted • • •). 

in order to have a non-trivial P. 

2.2.2 Reasons for instability when (3 > 2.2 

Novi^, for high values of P (i.e. (3 > 2.2) we faced difficulties in finding a solution. 
We believe that this has the following explanation: as P grows, the area of values 
where the potential U is negative, increases as well. 

If the parameter a is big, then the potential becomes a strong factor in reduc- 
ing the energy, as it is highly negative. Thus, for big a it is energetically favorable 
to decrease further the potential value. The latter happens by converting P func- 
tion in such a way, so as to be, as much of it as possible, inside the energetically 
favorable area (it is denoted by the horizontal lines in the figslHMl)- Thus, we end 
up to a final configuration with potential U « and virial relations can not be 
satisfied (see for example eqJT2ll as we have a sum of positive terms. 

On the other hand, for lower a the potential is no longer a strong factor for 
reducing the energy. In this case, it is energetically favorable to reduce the value 
of the {dpPy and {dpA^^Y terms as the system can gain more from these. The 
consequence is that P leaves the area of stability as its peak lowers in order to 
reduce the two terms above and there is only one possibility: to end up to zero 
energy, that is to say, the trivial configuration. 

The difference in the stable solutions we have found above, is that the values 
of P are such, that the potential is negative but not strongly negative while the 
changes on the terms {dpPy and {dpA^p)'^ can deform the field P in such a way, 
so that it can still be inside the favorable area. Then, the potential term has the 
possibility to change in such a way, so it can satisfy virial as well. 

2.2.3 Reasons for instability when (3 < (3crit 

In the following explanation we will use figlS & El In figlS one can observe that 
as we get closer to the critical value (3 = 3/-\/2 = Pcru, the P field tends to 
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Figure 6: The P+ and P- solution of eq.^, for which the potential becomes zero. 
Between these lines the potential gets negative values (see also hgE])- The plot is P vs. 

/3 for /? > fdcrit. = 



acquire everywhere the value P = Pq (location of the non-trivial minimum of the 
potential). This leads to greater values of dpP. 

The above has a reasonable explanation which can be found in figlGl As we 
get closer to Pcru, the space within the lines, where the potential can get negative 
values, becomes smaller. Under Pcru the potential can be either positive or zero 
(the latter for P = only). 

Observe the energy functional to be minimized: 



E = 27i pdp 







7^,id,A^r + {d,Pf + ^{eA^ + Mf + U{P) 
zp p 



(13) 



The main target of minimization is to "fix" all the above terms in order to have 
the minimum possible value for the energy. Below, we analyze the possible cases. 

In that case, all terms except for the potential term, can be either positive 
or zero. The potential term (as we saw in figlH]) can become negative for 
a range of values of P. But, as [3 decreases, this range becomes narrower 
(as we see in figE] as well as in figslSHU where this range is represented by 
the space between the two horizontal lines) and (dpP)^ increases. This 
happens because the energy functional tends to decrease its value through 
the negative values of the potential term. We believe that this can not 
continue for (3 very close to Pent due to the fact that the range we described 
above, becomes so small, that P tends to get everywhere a constant value 
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Figure 7: The plot exhibits the relation of the parameter a with respect to /5. We 
observe that as we approach the limit PcrU, we need an increasingly deeper minimum 
of the potential which is expressed through the fastly increasing value of a. The dashed 
vertical line signals the position of PcrU = 

(the value Pq, which makes the potential negative). But P must be zero 
at p = and p ^ oo, thus there will be a considerable increase in the 
[dpP)"^ term of the functional, and that makes the benefits of the negative 
value of the potential to go away, while the trivial solution P = becomes 
energetically favorable. 

• /5 < (^crit- In that case, the potential term can no longer become negative. 
It can be either positive or zero and because of the fact that the energy 
functional is a sum of five positive terms, it's reasonable to prefer the zero 
value which, at the same time, minimizes all the terms of the functional. 

From all the above, one can observe that the crucial difference between the 
above two cases of /?, is that in the {3 > jScrit case, the energy can have a minimum 
value (through the potential term) which corresponds to a non-trivial solution 
for P. The existence of negative values of the potential are the "way-through" 
that make this possible. 

3 The U{1)a X U{l)w model 

Stable strings of the above U{1) model could be useful in order to create vortex 
rings and study their stability. It could also be relatively helpful numerically, as 
we would have four fields for minimization in the energy functional instead of five 
we had in [TH|. This is not as easy as it might seem, since there are two instability 
modes. The vortex itself is not necessarily stable while forming a torus, and the 
latter has the tendency to shrink due to its tension. The first instability can be 
avoided in a U{1) x U{1) model as the one presented in [18], where a numerical 
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search for bosonic superconducting static vortex rings in a U{1)a x U{1)w model 
is done. There, the existence of straight strings is ensured for topological reasons. 
The superconductivity of the loop though, does not seem to prevent shrinking. 
The conclusion there is that current quenching takes place before stabilization. 
The model discussed in [18] is being described by the Lagrangian density 

^0 = -^^i - \w'^. + \D,^P\' + \D,<I>\' - UM, 1^1) (14) 

where the covariant derivatives are D^ip = d^ip + i^Af^ip, D^cp = d^cp + iqW^cj), 
the strength of the fields are F^i, = d^A^ — d^A^, Wfj,^ = dfjWu — d^W^, while e 
and q stand as the relevant U{1) charges. The potential U is 

^(101, 1^1) = - viY + f (i^r - viy + |i0n^p - (15) 

The vacuum |0| = fi 7^ 0, {ipl = 0, breaks U{l)w x U{1)a — ^ U{1)a, giving 
non-zero mass to W. The photon field stays massless. There, U{vi,0) = 0. The 
vacuum manifold A4 in this theory is a circle S'^ and the first homotopy group 
of M. is TTi{Ai) = 7ri(S'^) = Z which signals the existence of strings. In regions 
where |0| = 0, the field \ip\ is arranged to be non-vanishing and U{l)w x U{1)a 
U{1)]Y. Thus, U{1)a — ^ 1 and electric current fiows along regions with vanishing 
\(f)\. Hence, this theory has superconducting strings \19\. The vacuum of the 
theory leaves unbroken the electromagnetic U{1)a- For g^vf > g2V2 this vacuum 
is stable, while givf > (72^2 ensures that it is the global minimum of the potential. 
The mass spectrum is 

mA = 0, mw = qvi, ml = givl, ml = ^{gsv^ - g2vl) (16) 
3.1 The extended U{1)a x U{1)w model 

We intend to modify the above model by adding higher derivative terms of the 
fields and ip and find out whether such changes can stabilize the ring. By fol- 
lowing Derrick's scaling argument |2T], one can argue that terms such as l-D^^V'l^ 
or 1-0^01"^ or |-D^0p|D^^p could be helpful. Also, in an investigation of a similar 
model [201, the conclusions lead to the same path, in order to search for possi- 
bilities of stabilizing such solitons against radial shrinking. Thus, we have the 
following Lagrangian density: 

C = Co + c^\D^^\^ + c^lDf.ijl" + c^^\D^^\^\D^ij\^ (17) 

where c^, c^, c,/,^ constants. 

Configurations with torus-like shape, representing a piece of a U{l)w 1 
Nielsen-Olesen string, closed to form a loop, are of interest in this search. Thus, 
we will require to vanish on a circle of radius a (the torus radius) 0(p = a,z = 
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0) = 0. At infinity (p — oo, 2; — > 00), we have the vacuum of the theory. This 
translates to |0| — ^ vi^ \ip\ 0. The ansatz for the fields is: 



A(p,<^,2;) = — V 
W(p,(^,z) 



P 



Wp{p,z) p + W,{p, z) z 



where M, are the winding numbers of the relevant fields, p, 0, z are the 
cylindrical unit vectors and the phase 



9(p, z) = arctan ^ 



z 



p — a 



(18) 



We use cylindrical coordinates {t,p,ip,z), with space-time metric that has the 
form g^y=diag{l, —1, — p^, —1). We work in the AP = Q = W'^ gauge. We follow 
the ansatz of [I8|. A more general choice for would be 

0(p, (/^, z) = F{p, z)e'^mP,-)+ixiP,z) 

where x(p, 2;), an arbitrary function. But gauge invariance allows us to change 
— i> (f)e'^''^^\ where 6(x) an arbitrary space-dependent phase. We can choose 
fe(x) = — z) thus, gauge fixing removes the arbitrary function X- 
With the above ansatz, the energy functional takes the form: 



E 



2txvi I pdp 



dz 



2p 



+ {d,PY + {d-^Py + {d,FY + {d^FY + —{eA^ + NY + 

+ ({qWp + MdpQY + {qw, + Md.eY) + 

+ c^\^{dpFY + {d,FY+[{qWp + MdpQY + {'lWz + Md,QY)F^^ + 

+ c^\^[dpPY + {d^PY + ^{eA^ + NY^ + 

-f c^^l + {d^FY + ((gW^p + MdpQY + (gW^. + M9,e)2)F2j 



{dpPY + (5,P)2 + ^(eA^ + AT)^ 1 ;> + f/(F, P) 



(19) 
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with 

U{F, P) = - 1)' + j{P'- u'Y + ^-^F'P' - (20) 

where u = V2/V1. This is the energy functional we use for our computations. The 
conditions to be satisfied by the parameters become: 

91 > 92u'^ , 5-3 > 92U^ (21) 

The magnetic fields are 




dW, _ dWp 
dp dz 



while the currents associated with field, namely and jf and the total current 
X'^ out of these as well as the supercurrent X'^ associated with the field are 

I o p2 

I'^ - y Utr + Utr , X^^-^{eA^ + N) (22) 

where 

fp^-2qF^{qW, + Md,e), jf = -2qF\qW, + Md,Q) (23) 

Finally, in order to check our numerical results, we can derive virial relations. 
Below, we present the virial relations we use in our search. Consider the rescalings 
p p, z ^ Kz, ^ F, ^ P, A^^ A^, Wp^^^ K-Wp^^. By demanding 
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^ = when k—1 and if we define 



/■OO /-OO 



Ii — 2'KV\ I pdp / dz 



+ 



2P2p2 
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27rt'i / pdp I dz 



— (a,^^)^ + {d,Ff + {d,Pf + '-{eA^ + Nf + 



p2 
p2 



+ 



p4 



+ c<^v,<^ {dpPf{dpFf + ^(e^^ + Nf{dpFf + 



+ 



p2p2 

7 



+ 



(eA^ + Nf + F2(9,P)2 (giy^ + MdpQf + (gM^, + M^.e)^ ) + 



4 4 2 4 



we must have Ii + I2 = 0. We define the index V = and we want its 

value to be as small as possible. We can produce many other virial relations by 



assuming generally for a field 0, the "double" rescaling (j){x) K(j){jjix) and then 
demand 1^1^=1=^, = = |^|k=i=^. For example, we check our results through 
the following relations as well. Consider the following rescalings p — > p, ;s — > pz, 
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F« ^F,P,^ kP, A^^ ^ A^, Wp,,^ ^ Wp,,. We define 

1 



2'KVi 



poo /"C 

/ pdp I 

Jo J-^ 



dz 



+ iF" (^{qWp + MdpQf + {qW, + Md^Qf^ (Md,e{qW, + MdS)) 



+ 



+ 



+ cA 3{d,P)^ + 2{dpP)\d,Pf + '^{d,P)\eA^ + Nf) + 



+ 2F'' ( (a^P)' + (a,P)2 + ^{eA^ + TV)' I {Md,Q{qW, + M^.G)) + 
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oo 

dz 




oo 



— (a,^^)^ + {d,Pf + {d,Ff + ^(e^^ + Nf + 



p2 
p2 



+ 



+ 



p4 



{qWp + MdpQf + {qW, + Md.Qf 



4 4 2 4 



together with 



/■oo poo 

I5 — Anvi / pdp / 

-/o -/-oo 



{d,P)' + {d,Pf + —ieA^ + Ny + 



2P 



+ 2{d,Pf + 2{d,Pf + ^(eA^ + Nf + 4(a,P)^(a,P)^ + 



+ 



4p2 



((a,p)2 + (a.p)2)(eA^ + iv)4 + 



+ c^^|(a,p)2(a,p)2 + {d,p)\d,Ff + (a,p)2(a,p)2 + (a.pj^ia.p)^ + 



+ 



^{eA^ + Nf{{d,Pf + {d,Pf) 



POO P 

L 


00 

dz 




00 



n2 d2 
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where, as above, we must have I3 + = = + Iq. 

In the energy functional (fTOll . the terms that come from the |-Di0|^ extra term, 
are multiplied with c^. In fact, these terms are proportional to dpF and dzF. 
Thus, if one chooses (c0,c^,C0^) = (1,0,0), then there are more F-derivative 
terms in the functional. Energy minimization lowers these terms, something 
which one expects to lead to a thicker string. This is a wanted feature in order to 
stabilize the ring. This would have another consequence. The extra "F-terms", 
enforce the F field to stay away from its vacuum expectation value within a larger 
area. This means that F ^ inside a larger area. There, the potential becomes 



and its value increases because of the term gi{F'^ — 1)^/4, as F ^ 0. Minimiza- 
tion tends to make P u, which tries to compensate for that increase. The 
latter means that P will increase and this is another wanted feature since the 
supercurrent oc P^. 

On the other hand, in (fTOl l. the terms that come from the l-DjV'l^ extra term, 
are multiplied with c^. In fact, these terms are proportional to dpP and dzP- 
Thus, if one chooses {C(i,,c^,c^^) = (0,1,0), then the "derivative" terms of P 
become more. The minimization of them is expected to decrease the charge 
condensate P. This decrease is an unwanted feature. 

Finally, the terms that come from the \Di(j)\'^\Diip\'^ term, are multiplied with 
c^^. The consequences of the addition of this term can be seen if we observe that 
the extra terms are of the form 



where the dots represent the rest of the terms which are positive as well. We 
expect a stronger decrease of the charge condensate P. This is because dP < 1 
which means that (dP)"^ > (dP)'^. Thus, the need for minimizing the derivative 
terms of P becomes stronger than in the case of (c0,c^,C0^) = (0, 1,0). Apart 
from this, it is also the fact that, since dF ~ 1, the weight of the P-derivative 
terms is now greater than unity and this is another factor which would tend to 
make P — *• or, at least, smaller than in the case (c^, c^, c^^) = (0, 1, 0). 

Theoretically, the fact that the terms l-Dj^'l'^ and |P'i0p|Dj?/'p tend to shrink 
the charge condensate P, can also be seen from the virial relation Is + Jg = 
above. The integral Iq is negative (because P < u) and the addition of extra 
terms leaves it unchanged. On the other hand, is a sum of positive terms and 
the above two extra derivative terms rise the value of I^. Thus, the only way to 
satisfy that virial relation is either to increase g2 and/or to decrease P. If /s is big 
enough, then P will be enforced to become zero in order to satisfy the h + I^ = 
relation. 




(24) 



{dPf {dpFf + {d^Ff + ■ ■ ■ 



(25) 
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Figure 8: A typical plot of the initial guess we use, for the lowest winding state 
M = 1, iV = 1 on the z = plane. 

3.2 The extended U{1)a x U{1)w model: Numerical results 

We use the same minimization algorithm, as in [1^, to minimize the energy 
functional (flQl ). A 90x20 grid for every of the five functions is used, that is, 90 
points on p-axis and 20 on z. We begin with fixed torus radius a. Then, the 
configuration with minimum energy for this a, is found. Other values of a are 
chosen as well and the same process goes on until we plot the energy vs. the 
torus radius E{a). It would be very interesting to find a non-trivial minimum of 
the energy (on ttmin 7^ 0), which would correspond to stable toroidal defects with 
radius Omm- One crucial check of our results is done through virial relations. 
The initial guess (figl8|) we use for our computation is: 



tanh((p-a)2 + ^2)Af/2 
tanh(p^)(l - tanh((p - + z^) 




) 



2 



Mzcos'^ef {p-af + z^ 



q{p - a)2 \{p- a)2 + + (a2/4) 



W,{p,z) 



Mcos^ef (p-a)2 + z2 



) 




This initial guess also satisfies the appropriate asymptotics 
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Figure 9: A typical graph which exhibits the effects of the higher derivative terms on 
the scalar fields. The plot is for the z = plane. Solid line is for {c^,c^,c^^)={0,0,0), 
dotted line for (c,^, c^, C(^^)=(0, 1, 0), dashed for (c^, c^, C(^^)=(l, 0, 0) while dashed and 
dotted for (c,^, c^, C0^)=(1, 1, 0). In the case (c^, c^, C0^)=(0, 0, 1), P is trivial. Param- 
eters in this figure are (51, 92, 53, e, g, u, f i, M, A^)=(14, 12, 14, 6, 2, 1, 7.5 • 10"^, 1,1). 

• near p = 0: 

F^O, P~p^, A^^p'fiz) (26) 

• near {p = a, z = 0): 

Fr^p^'/"^, Wp = = W, (27) 

• at infinity: 

F ~ l-0(e-^), P~C(e-v/^) 

W, i^5,eU + 0(e-v^), ly, a,GU + 0(e-v^) (28) 

q q 

where p = {p — a)^ + z"^. 

Based on [18j, we search in parameter areas where e acquires relatively large 
values, but they are interesting as it concerns the possible stability of the loop. 
The reason was analyzed in that paper and stems from the need to have string 
thickness greater than the penetration depth as well as strong supercurrent. 

The results confirm our expectations stated previously. For example, in figlHl 
one can observe that when the extra term is then the charge condensate 

P decreases (dotted line in figl9|). In this case, the consequence is the reduction of 
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Figure 10: The effect of the higher derivative terms on the energy and the supercurrent 
when these are plotted vs the radius of the torus. Sohd hne is for (c^, c^, C0^)=(O, 0, 0), 
dotted line for (c^, c^, C(^^)=(0, 1, 0), dashed for (c,^, c^, c<^^)=(l, 0, 0) while dashed 
and dotted for {c^,c^,c^^)={l,l,0). For (c<^, c^, C0^)=(O, 0, 1), the supercurrent 
as well as the charge condensate, are trivial. Parameters in this figure are 
(51, 92, 53, e, q, u, v^M, N)={U, 12.5, 14, 10, 2, 1, 7.5 • lO'^, 1, 1). 



the supercurrent (dotted line in fig JTOl) when compared to the case of the original 
model without extra terms ((c<^, c^, c,^^) = (0,0,0), see solid line in fig JTOl) . On 
the other hand, when the extra term is \Di(j)\^, then F widens and this leads to 
the broadening of P as well. The latter increases (compare dashed and solid lines 
in figl9|) and the supercurrent increases too (dashed line in figJTOl). Finally, when 
the extra term is then P = for the values of g2 we use. In general, 

the charge condensate can be non-trivial for higher g2. This happens due to the 
term g2{P^ — m^)^/4 of the potential. When g2 grows, P tends to reach u. This 
is also numerically observed. 

One can combine two extra terms to see what happens. In example, we add 
both \Di4>\^ and \Diip\'^, (case (c,^, c^, C(^^) = (l, 1, 0) in figl9|). This results to the 
addition of the "favorable" F-derivative terms, but also P-derivative terms would 
be present. This translates to the growth of P but not as much as in the case 
(c0, c^, C(/,^)=(l, 0, 0). We also observed that the combination of either l^j^l^ or 
lAVI^ or 6ven both, with the \Di(f)\'^\Di'ip\'^ term, leads to shrinking of P because 
of the strong action of the last term. After the theoretical and numerical analysis, 
we conclude that the most "interesting" extra term is |P'i0|'^. 

It is clear that current quenching is present here as well (figHOj). The extra 
term lA^I^, can increase the supercurrent and can make the penetration of the 
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magnetic field more difficult, as the string increases its diameter, but this increase 
is not enough in order for the ring to stabilize. For the shake of research, we also 
tried higher values of c,^ in order to make the favorable term more significant. We 
also tried higher values of e, but the ring could not be stabilized in a non-trivial 
radius. 

3.3 Discussion 

The most crucial terms of the energy functional which could provide for the 
stability of the ring in a non-zero radius a, are: 




B = / pdp dz —{eA^ + Nf 

Jo J-oo 

These two terms have an explicit total 1/p behavior which helps them to increase 
as the torus radius decreases. The problem is that they are not increasing at a 
satisfactory rate in order to overcome all the rest terms of the energy which are 
decreasing with p. This is also numerically observed. Under some radius a, the 
ideal would be to have a strongly increasing charge condensate P. Then, as the 
radius decreases and at the same time P increases, the above two terms would 
start to increase with sufficient rate in order to lift the energy of the system. The 
magnetic term would increase because as P'^/p^ increases, \A^p\ \A^^^J —>■ 
N/e. The latter would make the derivatives of A^ (and Ba as well) to increase 
as p decreases. 

This is what we tried to do here, especially with the help of the |-Di0|^ extra 
term. The charge condensate became more robust but that was not enough. This 
supports the conclusion of \i8\ which states that very high values of supercurrent 
are needed for stabilization. It seems that such highly increasing currents can 
not be produced, despite the help of extra terms. In fact, numerical details reveal 
that the rate of increase of the terms A, B above, is A{A + B) ~ 10~^, while the 
rate of decrease of all the rest terms is A{E — {A + B)) ~ 3 ■ 10~^. This means 
that the rate of increase of A, B should be ~ 30 times bigger. That case would 
require P >u which is something that can not satisfy Is + Jg = virial relation. 
But even if that was possible, current quenching would be another "obstacle". 

4 Conclusions 

We studied a U{1)a model with a GL potential with a cubic term added to it. 
After the numerical analysis we did, we came to the conclusion that for /3 < /^crit 
we find no non-trivial solution. For 2.13 </? < 2.20 we get non-trivial solutions 
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which have the profiles we present in figsjSHU Over [3 = 2.20 we have no non- 
trivial solutions. We analyse and explain our results in these cases. 

The form of the potential we have in this search, can be found in condensed 
matter physics as well. On the other hand, one could try to make a loop out of 
the straight string studied above. We did this, but it was difficult to study mainly 
due to the fact that there are two instability modes, one having to do with the 
defect itself and another which has to do with the loop that tends to shrink due 
to its tension. The former instability is excluded in [18] for topological reasons. 

Thus, we are based on the model of pE| and analyze an extended version of 
it, by adding higher derivative terms in order to check whether they can stabilize 
the superconducting ring or not. Although the l-Dj^l"^ term is helpful on that 
direction, it turns out to be insufficient and current quenching prevails. Finally, 
we discuss what one would need for a stable ring. This discussion in combination 
with the results of [I8], seems to exclude the possibility of existence of such vortex 
rings in this model. 
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